###***Create lists (data storage for each bin)***###

r = []              #radius of box
deltaA = []         #Area of box
m = []              #mass in box
sigma = []          #surface mass density
R = []              #Effective r**2 factor for moment of inertia
I = []              #moment of inertia of box
w = []              #orbital speed of box
R_M = []            #Radius in units of Mars Radii
R_S = []            #Radius in units of Saturn Radii
Ncol = []           #Number of collisions with Mars in timestep
Ecc = 0.09341233    #Eccentricity at Mars of Heliocentic disk
Inc = 1.85061       #Inclination of Heliocentric disk at Mars
ecc = 0.02#[]       #Eccentricity of disk orbiting Mars as a function of time
inc = 0.01#[]       #Inclination of disk orbiting Mars in time
t = []              #time
Mars_M = []         #Mass of Mars in time from datafile
Mars_R = []         #Radius of Mars in time
Mars_DeltaM = []    #change in mass of mars in time
NImp = []
coll = []
spread = []
vesc = []           #Escape speed b/w Mars and impactor
vsw = []            #Swarm speed (relative speed) b/w "
v = []              #Collision speed b/w "
probability = []    #Distribution for probability a box will be hit
nu = []             #Viscosity value
signew = []         #New surface mass density
u = []              #Used in visc_spread routine
delm = []           #Used in Disk_Restructure
w_new = []          #Used in Disk_Restructure
Mass_to_Mars = []   #Mass accreted to Mars for each timestep n
DiskSum = []
MarsAcc = []
Sat_M = []
Sat_R = []
m_new = []
moons = []
vec = []
Sat_M = []
Sat_Rad = []
Sat_r = []
Sat_r_default = []
Sat_v = []
Sat_w = []
Sat_e = []
X = []
S = []
rhill = []
Torque_to_disk = []
Sat_tstep = []
Laplace_Coefficient = []
DLaplace_Coefficient = []
